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ABSTRACT 

We combine Hubble Space Telescope images from the Great Observatories Origins Deep Survey with 
archival Very Large Telescope and Keck spectra of a sample of 11 X-ray selected broad-line active 
galactic nuclei in the redshift range 1 < z < 2 to study the black hole mass - stellar mass relation 
out to a lookback time of 10 Gyrs. Stellar masses of the spheroidal component (Msph,^) are derived 
from multi- filter surface photometry. Black hole masses (Mbh) are estimated from the width of the 
broad Mgll emission line and the 3000A nuclear luminosity. Comparing with a uniformly measured 
local sample and taking into account selection effects, we find evolution in the form MBH/Msph,^cx 
(1+z)^"^^^^"^^, in agreement with our earlier studies based on spheroid luminosity. However, this result 
is more accurate because it does not require a correction for luminosity evolution and therefore avoids 
the related and dominant systematic uncertainty. We also measure total stellar masses (Mhost,*)- 
Combining our sample with data from the literature, we find MBH/^host,^cx (1 + z)^-^^^^-^^, consistent 
with the hypothesis that black holes (in the range Mbh lO^^^M©) predate the formation of their 
host galaxies. Roughly one third of our objects reside in spiral galaxies; none of the host galaxies reveal 
signs of interaction or major merger activity. Combined with the slower evolution in host stellar masses 
compared to spheroid stellar masses, our results indicate that secular evolution or minor mergers play 
a non-negligible role in growing both BHs and spheroids. 

Subject headings: accretion, accretion disks — black hole physics — galaxies: active — galaxies: 
evolution — quasars: general 



1. INTRODUCTION 

Active Galactic Nuclei (AGNs) are thought to repre- 
sent an integral phase in the formation and evolution 
of galaxies during which the central supermassive 
black hole (BH) is growing through accretion. The 
empirical relations between BH m ass (Mbh) and the 
prope r ties of the host galaxy f e.g.. iFerrarese fc MerrittI 
mm IGebhardt et al.l I2OOOI : IMarconi fc HuntI 120031 : 
IHaring fc Rixl l2004f ) have been explai ned by a com- 
binat i on of AGN f eedback fe-p^.. [Volonteri et al.l 
[2OOI ICiotti fc Ostrikerl l200'7l : iDi Matteo et al.l 120081 : 
iHopkins et al.l l2009f ) and hierarchical assembly of 
Mbh and stellar mass throu gs:h galaxy merging (e.g., 
[Pen^[2007l : Ijahnke et al.ll20TTh . 

The great interest in the origin of the scaling 
relations is reflected in the flood of observational 
studies, f ocusing on their cosmic evolution (e.g., 
Treu et al.l [20 04: Wal ter et al. 2004; Shiel ds et al l 
McLure et al. 2006; Peng et al.. .2006: .Woo et al.l 



2006; 



Salviander et al. 2007: iTreu et al .11200 TtTWoo^^alJ 



2006; 



2008; 



Jahnke et al. 2009; Bennert et al.J 120101 : iDeca rli et al 
2010; Merloni et al. 2010), with the majority pointing to 
a scenario in which BH growth precedes bulge assembly. 

However, many high-redshift studies to date are based 
on monochromatic Hubble Space Telescope (HST) imag- 
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ing, determining only the luminosity of the host spheroid 
and not its stellar mass. This is acceptable at z ~ 0.5 
(e.g.. iBennert et al.ll2010f ). where the stellar populations 
of bulges are fairly well known and their luminosities can 
be passively evolved to zero redshift with uncertainties 
smaller than other sources of error. In contrast, at 2; > 1, 
the stellar populations of bulges are an uncharted terri- 
tory, particularly for AGN hosts which are believed to be 
connected with major mergers and ma y have undergone 
recent episodes of st ar formation (e.g.. iKauffmann et al.l 
120031 : iSanchez et al.h200 4). The uncertainty on the con- 
version from observed luminosity to equivalent 2: = 
luminosity can b e comparable to the evolutionary signal 
(e.g., Peng et alJE006). ^ 

An exception is the study by iMerloni et al.l ()2010l ) 
who estimate total stellar masses (Mhost,^) and AGN lu- 
minosities by fitting spectral-energy distribution (SED) 
models to multi-band data from the rest-frame ultra- 
violet to the rest-frame mid-infrared for a sample of 
89 broad-line AGN (BLAGN) hosts at 1 < z < 2.2. 
Estimating Mbh from broad MgH emission, they find 
that black holes of a given mass reside in less mas- 
sive hosts at higher redshi ft with a modest evolution- 
ary slope. However, Merlon Tet al.l (2010) are unable to 
distinguish between Mhost,^ and the stellar mass of the 
central bulge component of the host (Msph,^). Such a 
difference may be important when studying the evolu- 
tion of the scaling relations: there are indications that 
the relations between Mbh and total host-gal axy luminos- 
ity (Benne rt et al.ll2010t) and stellar mass (Jahnke e t al.l 
1^09) may not be evolving, or at least not as rapidly as 
the relations between Mbh and spheroid properties. 

In this paper, we study the cosmic evolution of the 
MBH-Msph,^ and MBH-Mhost,^ relations for a sample of 
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11 BLAGNs {1 < z < 2; lookback time: 8-10 Gyrs) 
selected from the Great Observatories Origins Deep Sur- 
vey (GOODS) fields, taking into account selection effects. 
Msph,^ and Mhost,^ are derived from the deep mult i- filter 
HST images. Mbh is estimated using the width of the 
broad Mgll emission line, measured from existing spectra 
and the 3000A nuclear luminosity. We use a local com- 
parison sample of Seyfert-1 galaxies (jBennert et al.|[2QTTf ) 
for which all relevant quantities were derived following 
the same procedures adopted for the distant sample to 
minimize potential systematic bias. Our strategy allows 
us to address two major hmitations of previous studies: 
eliminate uncertainties due to luminosity evolution and 
determine the evolution of the spheroidal component of 
the host. 

Throughout the paper, we assume a Hubble constant 
of i^o = 70kms-i Mpc-\ Qa = 0.7 and = 0.3. Note 
that all magnitudes are AB. 

2. DATA ANALYSIS 
2.1. Sample Selection 

The hi^h-redshift sani ple consists of AG Ns in GOODS- 
N (Treister et al.' ^2009) and GOODS-S (jTrouille et al.l 
I2OO8; Silverman et al. 2010), selected based on their 
X-ray emission using the Chandra Deep Field North 
and South (CDF-N/CDF-S) survey and spectroscopi- 
cally confirmed to be BLAGNs. We select all 11 objects 
within 1 < z < 2 for which archival Very Large Tele- 
scope (VLT) and Keck spectra covering the broad Mgll 
line exist (Table [I]). 

By design, all objects have deep HST/ Advanced 
Camera for Surveys (ACS) images in four different 
broad-band filters (B=F435W, V= F606W, z=F775W, 
z=F850LP) (jGiavalisco et al.ll2004f ). Color images are 
shown in Figure [TJ The total exposure times range be- 
tween 5,000 and 25,000 sec, depending on the filter and 
the image region. The reduced data are taken from the 
v2.0 data release.^ The spatial resolution is approxi- 
mately O'.'l full- width- at-half- maximum (FWHM), which 
at z = 1.3 (our average redshift) corresponds to 0.84 kpc; 
thus, our data have higher spatial resolution than Sloan 
Digital Sky Survey (SDSS) images at z = 0.05 {VfA = 
1.37 kpc). Overall, the AGN host galaxies look like typ- 
ical ellipticals or spirals, without any signs of merger ac- 
tivity. 

Our local comparison sample consists of 25 Seyfert- 
1 galaxies selected from SDSS (0.02 < z < 0.1; Mbh> 
IO^Mq) for which all relevant quantities were derived fol- 
lowing the same procedures adopted for t he distant sam- 
ple to minimize potential systematic bias (jBennert et al.l 

[mil). 

2.2. Surface Photometry 

We perform two-dimensional surface photometry using 
the code "Surface Photometry and Structural Modeling 
of Imaging Data" (SPASMOID) developed by one of us 
(M.W.A.). The code allows a joint multi-band analysis 
of surface brightness models, thu s superceding the func- 
tionality o f GALFIT (iPeng; et al1 [2002). and is described 
in detail in lBennert et al.l ()2011 f). The point-spread func- 
tion (PSF) of the HST/ ACS optics is modeled using the 

^ http://archive.stsci.edu/pub/hlsp/goods/v2/ 



closest bright star to a given object. We impose a Gaus- 
sian prior on the AGN colors with the mean given by 
the qu asar composite spectrum from lVanden Berk et al.l 
(|2001f ) redshifted to the AGN redshift and with a cr of 
0.2 mag. We model the host galaxy by eithe r a sin- 
gle Ide VaucouleursI ()1948f ) profile or by a de Vau couleursI 
(jl948f ) profile plus an exponential profile to account for 
a disk. Depending on the images and residuals, we 
decide whether a given object is best fitted by three 
components (PSF, spheroid, disk) or two coni p onent s 
(PSF + spheroid^), as described by iTreu et al.l (I2007D : 
Benne rt et al.l (j201Ql I2OIII ) . A disk component is evident 
in 4/11 objects. In all four cases, we can only clearly de- 
tect the bulge in the z band. 

To probe the reliability of our AGN-host-galaxy de- 
compositions when using the blue restframe wavelengths 
covered by the GOODS images, we tested the effect of 
bandpass shifting. Given the host-galaxy morphology 
and level of activity of the sample studied here, a lo- 
cal sample of Seyfert gala xies is a suitable co mparison 
sample for this test. (Sch awinski et all l2010l also con- 
cluded that "moderate luminosity AGN host galaxies at 
z 2 and z ~ are remarkably similar".) We thus 
repeated the analy sis of our local sample of AGN host 
galaxies (Bennert et al.l [2QTTI ). but now using only ug 
SDSS photometry (instead of griz). We are able to re- 
cover the photometry of the bulge and point source to 
within 0.1 mags, i.e. smaller than our adopted system- 
atic uncertainty, demonstrating that bandpass shifting is 
not a concern within our level of precision. Moreover, 
this is a conservative estimation, since the GOODS im- 
ages at z 1.3 not only cover wavelengths comparable 
to ug rest frame (F775W and F850LP), but additionally 
also shorter wavelengths (F606W and F435W), thus ef- 
fectively providing more information to disentangle point 
source and bulge. Furthermore, as already pointed out 
above, the GOODS images of 2: 1.3 objects have even 
higher resolution than SDSS images at 0.5. 

2.3. Stellar Mass 

From the resulting magnitudes (Table [2j), stellar 
masses are estimated using a Bayesian stellar-mass es- 
timation code (Auger et al. 2009) assuming a Chabrier 
initial mass function (IMF) (Table [1]). We impose con- 
servative uncertainties of 0.3 dex on the masses of the 
bulges (disks) for the bulge-dominated (disk-dominated) 
hosts. The masses for the bulge components of the disk- 
dominated hosts are estimated by using the z-band mass- 
to-light ratios of the bulge-dominated hosts in our sam- 
ple that are at similar redshift s; we therefore add in 
quadrature a 0.3 dex uncertainty, yielding a total stel- 
lar mass uncertai nty of 0.4 dex for these objects. For 
two of our objects. [Schawinski et al.l ()2010f ) report stellar 
masses based upon template fits to the integrated light. 
Our results a gree within the un certainties (assumed to 
be 0.2 dex for I Schawinski et al.l[20To[ ). 

2.4. BH Mass 

Black hole masses are estimated via the empirically 
calibrated photo- ionizati on method (^Virial method") 
(e.g., iWand el et al. 119991 : IVestergaard PetersonI 12006^ 
iMcGih et al. 2008), by combining the FWHM of the 
broad MgIIA2798A emission line and the 3000A AGN 
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Fig. 1. Deep HST/ACS color images (B, V, i, z), 3'' xS" 
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TABLE 1 
BH Masses, 



AND Stellar Masses 



ID 


R.A. 


Decl. z 


Ref. 


FWHMMgii 


AL3000 


Mbh 




^disk,* 




(J2000) 


(J2000) 




(kms-l) 


(10^4 ergs-1) 






(1) 


(2) 


(3) (4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 



J033252-275119 
J033243-274914 
J033239-274601 
J033226-274035 
J033225-274218 
J033210-274414 
J033200-274319 
J033229-274529 
J123553+621037 
J123618+621115 
J123707+622147 



ID. 88 
ID. 24 
ID. 09 
ID. 50 
ID. 17 
ID. 91 
ID. 36 
ID. 98 
ID. 13 
ID. 58 
ID. 46 



ID.8 
ID.2 
ID.8 
ID. 5 
ID.8 
ID. 9 
ID. 7 
ID. 9 
ID. 3 
ID.O 
ID. 9 



1.227 
1.900 
1.220 
1.031 
1.617 
1.615 
1.037 
1.218 
1.371 
1.021 
1.450 



S10/S04 
T09/S04 
T09/S04 
S10/S04 
S10/S04 
S10/S04 
S10/L05 
S10/M05 
T08/W04 
T08/W04 
T08/W04 



16208 
16381 
4344 
2430 
4744 
5852 
3602 
5308 
5441 



10654 



0.43 
1.77 
5.38 
9.51 
1.64 
2.02 
1.08 
4.33 
2.32 
1.19 
1.57 



9.01 
9.31 
8.39 
8.00 
8.22 
8.45 
7.90 
8.52 
8.41 
8.49 
8.91 



9.83 
10.64 
10.54 
9.53 
10.61 
10.45 
9.62 
10.71 
9.99 
9.29 
10.74 



10 



10 



50 



75 



84 
95 



Note. — Column 1: target ID (based on R.A. and Decl.) . Column 2: Righ t Ascension. Column 3: Declination. Column 
4: redshift (taken from Team Keck Redshift Survey (TKRS) CWirth et al."2004') for GOODS-N). Column 5: reference for cat- 
alog f rom which objects were selected/ reference for origin of spectra. S10= Silverman et al. (2010), T09=Treister e t al.l 
(|2QQ9|) T08=Trouille et al. C2Q08); S04=Szoko ly_et al, (2004) , L05 =Le Fevre et al. (2005), MOS^MignoH et al. (2005^ 
WQ4= rwTr'th et al.l (120041) . Column 6: FWHM of broad Mgll. Column 7: rest-frame luminosity at SOOOA (fiducial error 
0.1 dex). Column 8: log Mbh /Mq (uncertainty: 0.5 dex). Column 9: stellar spheroid mass log M^p^^^^ /Mq (Chabrier 
IMF; uncertainty 0.3 dex for ellipticals, 0.4 dex for spirals). Column 10: stellar disk mass log Mdisk,*/-^© (Chabrier IMF; 
uncertainty 0.3 dex). 



TABLE 2 
Results from Surface Photometry 



Object 



PSF 



Spheroid 



Disk 





B 


V 


i 


z 


B 


V 


i 


z 


B 


V 


i 


z 




(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


J033252-275119 


25.03 


23.84 


23.40 


23.47 








23.86 


24.12 


23.49 


22.84 


22.18 


J033243-274914 


22.53 


23.09 


22.98 


23.10 


25.30 


24.68 


23.67 


23.37 










J033239-274601 


21.33 


21.08 


21.14 


21.38 


24.17 


24.32 


23.06 


22.10 










J033226-274035 


20.26 


20.04 


20.04 


20.01 








23.27 


22.22 


21.87 


21.22 


20.63 


J033225-274218 


25.38 


23.87 


22.84 


22.50 


25.08 


24.61 


23.84 


23.08 










J033210-274414 


24.24 


23.19 


22.61 


22.74 


24.13 


24.10 


23.55 


22.96 










J033200-274319 


22.81 


22.42 


22.45 


22.47 


43.90 


25.33 


24.00 


23.04 










J033229-274529 


21.24 


21.31 


21.42 


21.76 


23.80 


23.32 


22.49 


21.66 










J123553+621037 


22.55 


22.12 


22.15 


22.33 








24.04 


23.62 


23.20 


22.45 


21.83 


J123618+621115 


22.77 


22.27 


22.62 


22.81 








23.87 


23.15 


22.37 


21.40 


20.79 


J123707+622147 


23.16 


22.97 


22.69 


22.55 


24.95 


24.66 


23.46 


22.57 










Note. — Column 


1: target ID. Columns 2-5: extinction-corrected B, V, i, and 


z PSF magnitudes 


(uncertainty: 0.2 mag). 



Columns 6-9: extinction-corrected B, V, i, and z spheroid magnitudes (uncertainty: 0.2 mag). Columns 10-13: extinction-corrected 
B, V, i, and z disk magnitudes (uncertainty: 0.2 mag). 



+0.47 log 



3000 



10^4 ergs-1 



The AGN luminosity is derived from the PSF magnitudes 
in the filter closest to rest frame 3000A, and extrapolated 
based o n the assumed AGN SED of iVanden Berk et al.i 
([2001) (E2 Table [J). 



The nominal uncertainty of Mbh using this method 
is 0.4 dex. However, for some spectra, the low signal- 
to- noise (S/N) makes the FWHM measurements uncer- 
tain by up to ~50%, conservatively estimated. Moreover, 
the spectra are not of sufficient quality to remove the Fe 
emission which can result in overestimating the width of 
MgH by up t o 0.03 dex (in FWH M, McGi niral][2QQ8l : 
see, however. iMerloni et al.|[2QTQ[ ). We therefore adopt 
an uncertainty of 0.5 dex. Note that while we used uni- 
form priors for both the 3000A luminosity and the black 
hole mass in our analysis, employing more informed pri- 
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TABLE 3 

Evolving Mbh - ^* Scaling Relations 



Model 


a 


/3 


7 


a 


spn,* 


1.09 


1.96 ±0.55 


-0.48 


0.36 ±0.1 




1.12 


1.68 ±0.53 


-0.68 


0.35 ±0.1 


Kost,* 


1.12 


1.15 ±0.15 


-0.68 


0.16 ±0.06 




1.12 


1.11 ±0.16 


-0.68 


0.17 ±0.07 


Note. — '^Fitted only using the 11 objects pre- 
sented here. ^Fitted using the objects presented here 
and the objects from 'Merloni et al.| (|201Q|). '^Fitted 
only using data from Merloni et al.| (|2Q1Q[). 



ors from the quasar luminosity function of lRichards et al. 
(2006) or the black hole mass function of et al.l 

(|2010i) yield negligible changes to our inference. 

2.5. Mbu-M^ Evolution 

Foll owing and expandin g on work bv lTreu et al.l ([2007') 
and Bennert et al.l (|2010f ). we model the evolution of the 
offset of the MBH-^sph,^ and MBH-^host,^ scaling rela- 
tions by assuming a model of the form 

logMBH -S = a [logM* - 10] + f3\og [1 + z] + 7 + a 

where a is the slope of the relations at z = and is 
assumed not to evolve, 7 is the intercept of the rela- 
tions at z = 0, and a is the intrinsic scatter which is 
also assumed to be non-evolving. Here /3 describes the 
evolution of the scaling relation (with /3 = implying 
no evolution). We im pose J- function priors of a = 1.09 
(|Bennert et aIll201lD and 7 -Q-fS for the Mgph,* rela- 
tion and a = 1.12 (jHaring fc Rixll2004f ) and 7 ='-0.68 
for the Mnost,* relation; the priors on 7 were determined 
by fitting to the local AGNs from Bennert et al. (2011) 
while keeping the slope fixed to the noted values. A 
normal distribution prior is used for the intrinsic scat- 
ter with mean 0.4 and variance 0.01 and we employ a 
broad uniform prior for p. (Note that, strictly speak- 
ing, the variable a accounts for both the intrinsic scatter 
in the relationship and the (much smaller) uncertainty 
on 7.) We us e the z = 1-0 — 1.2 'elliptical' stellar mass 
function from lllbert et al.l (|2010[ ) to place priors on the 
stellar masses. Furthermore, we include a prior on the 
black hole masses that models our selection effects by 
using a hard cutoff at the low mass end. This cutoff is 
determined from the data and models lower limit of black 
hole masses observable in each considered set of data. 

The relation above is first fitted using the 11 galax- 
ies in this sample. The lower limit for the black hole 
masses assumed for t he high redshift objects is lO^'^M©. 
iMerloni et al.l (|2010f ) have independently tried to infer 
the evolution of the MBH-^host,^ relation, but their anal- 
ysis is somewhat different than ours (e.g. IMFs, local 
comparison samples, definition of offset, treatment of up- 
per limits and sele ction effects). We ther efore also fit the 
relation using the IMerloni et al.l ()2010l ) data (adjusted 
to a Chabrier IMF), and we impose a limiting black hole 
mass of 10^"^ for these data. The results of our inference 
are shown in Table [H Given that the different fits to the 
MBH-^host,^ relation (Merloni et al. data only, our data 
only, both combined) result in the same P within the un- 
certainties, we adopt the one for the combined sample in 
the following. 



3. RESULTS AND DISCUSSION 

4/11 AGNs are clearly hosted by late- type spiral galax- 
ies, while the rest seem to be spheroid dominated. Keep- 
ing in mind the small- number statistics, the fraction 
of disk-dominated host galaxi es (36 zb 17%) is lower 
than what has been found by iSchawinski et al.] ()2010l ) 
(80 ± 10%) for 20 X-ray selected AGNs at a comparable 
redshift (1.5 < z < 3) imaged by HST/Wide Field Cam- 
era 3 (WFC3; F160W) with 1-2 orbits integration time. 
One difference is that our objects have higher X-ray lu- 
minosities (0.5-8keV; 43.5 < logL^ < 44.5, mean=44.2, 
compared to 42 < \ogLx < 44, mean=43.1) which might 
explain why we find a larger fraction of elliptical host 
galaxies. 

Interestingly, none of the objects shows clear signs 
of interactions or merger activity, while at redshifts of 
z = 0.4 — 0.6, 32±9% of Se yfert-ls are host ed by interact- 
ing/merging galaxies ( Bennert et al.|[2Q10f ). However, we 
cannot exclude that some of these low sur face-brig:htness 
features might have been missed (see, e.g. Bennert e t al.l 

2008) . Schawinski et al. (2010) also do not report inter- 
actions/mergers but their images are significantly shal- 
lower than ours. Star- forming galaxies at a redshift of 
z ~ 2, on the other hand, sho w a 33^6% fraction of 
interacting or merging systems (Forster-Schreiber e t al.l 

2009) . Schawinski et al. (2010) interpret their high frac- 
tion of spiral galaxies as a sign that secular evolution 
may play a non-negligible role in growing spheroids and 
black holes. Our findings, including the lack of merger 
activity, are consistent with such a scenario. 

Figure [2] (left) shows the MBH-^sph,^ relation, in- 
cluding a s ample of 18 in active galaxies and the local 
AGNs from iBennert et all (,2011). In Figure [2] (right), 
we show the MBH -Mhost,^ relatio n , aga in including the 
local AGNs from IBennert et al.l (l201lD and a ddition- 
ally, the 89 AGNs from IMerloni eFaH (|20Tol ) (10/89 
with upper limits only; subtracting 0.255 dex to con- 
vert their total stellar masses from Salpeter to Chabrier 
IMF, Bruzual Charlotll2003f ). Note that for ah com- 
parison samples, Mbh were estimated using the same 
recipe adopted here. 

In Figure m we show the offset in log Mbh as a function 
of constant Mgph,^ (left panel) and Mhos t,^ (r ight panel) 
with respect to the z = relations (see §2.5p . For com- 
parison, the offset in log Mbh as a function of constant 
stellar spheroid lumino sity (left panel) and t otal luminos- 
ity (right panel) from IBennert et al. l (20101) is overplot- 
ted. Taking into account selection effects ( §2.5p . we find 
significant evolution in MBH/Msph,^((x (1 + zY'^^^^-^^), 
consistent with (but with larger uncertainties) what we 
reported previously for the evolution of the Mbh — ^sph 
relation (MBn/i^sphOC (l + z)^ '^^^-^; Bennert et aIll201Qf ). 
The agreement between the stellar mass and luminosity 
evolution suggests that the passive luminosity correction 
is appropriate, although modeling luminosity evolution 
rather than stellar masses may increase t he scatter. 

For total stellar masses, including the IMerloni et al.l 
(|2010l ) data, the evolutionary trend can be described as 
MBH/Mhost,^oc (1 + zY'^^^^-^^ , in agreement with what 
has been found by ^ Merloni et al.| pOlO) within the un- 
certainties. This evolution is slower than the one for 
spheroi d masses (0 = 1-96 ± 0.55) in line with recent 
studies (|Jahnke et al.l 120091 : [Bennert et al.ll2010f ). It in- 
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dicates that the amount by which at least some of the 
distant AGN host galaxies have to grow their bulge com- 
ponent in order to fall on the local BH mass scaling 
relations is contained within the galaxy itself. It can 
thus be considered as another evidence that secular evo- 
lution and/or minor mergers play a non- negligible role in 
growing spheroids through a redistribution of stars from 
disk to bulge. The deduced evolution is either in line 
with or slightly faster than what has been predicted by 
theoretical studies ( for a detailed comparison, see, e.g. 
iBennert et aL.2Ql(l iLamastra et al. 2010). 

4. CONCLUSIONS 

We determine spheroid and total stellar masses for the 
host galaxies of 11 X-ray selected BLAGNs (1 < 2: < 2) 
in GOODS. In combination with Mbh estimated via 
the virial method from the broad Mgll emission line 
as measured from archival VLT and Keck spectra and 
the 3000A nuclear luminosity, we study the evolution of 
the Mbu-M^ scaling relation out to a lookback time of 
10 Gyrs. Using a uniformly measured local comparison 
sample and taking into account selection effects, we find 
evolution of the correlations consistent with BH growth 
preceding galaxy assembly, confir niing and extending 
the results of previous studies (e.g.. [ Merloni et al.|[2QTQ[ 
iDecarh et al.ll2010l : iBennert et al.ll2010[ ). 

Our results show that a significant fraction (4/11) of 
AGNs at z=l-2 are hosted by spiral galaxies. None of the 
galaxies show evidence for recent major merger interac- 
tion, contrary to the general assumption that BHs and 
spheroids grow predominantly through major mergers, a 
scenario which might hold true only for the most lumi- 
nous AGNs. The evolution we find for the Men-total 
stellar mass relation is slower than the on e for spheroid 
stellar masses in line with recent studies (|Jahnke et al.l 
l2QQ9l : IBennert et al.llMol ). Combined, our results indi- 
cate that secular evolution and/or minor mergers play a 
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non- negligible role in growing both BHs and spheroids. 

Our study demonstrates the feasibility of obtaining 
stellar masses of AGN host galaxies out to lookback 
times of 10 Gyrs based on deep multicolor HST photom- 
etry. This approach has the great advantage of being 
independent of the luminosity evolution correction - the 
dominant source of systematic uncertainty in previous 
studies at compara ble redshifts (e.g. .Peng et al.l I2QQ6I : 
IBennert et al.|[2QTQh . Furthermore, we can distinguish 
between spheroid and total host galaxy mass, whic h is 
not possible based on SED fitting (e.g., Merloni et al.l 
[2010). 

Sample size is a major limitation of this work, allowing 
us to constrain only average evolution and preventing us 
from investigating, e.g, mass-dependent trends or corre- 
lations between evolution and morphology. Follow-up of 
BLAGN hosts imaged by existing and upcoming multi- 
color HST surveys (e.g. CANDLES) is needed to gather 
larger samples and address theses remaining issues. 
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Fig. 2. — Left panel: MeH-^sph,* relation for our sample (red pentagons; green circles: if fitted by spheroid plus disk), local BLAGNs 
(black circles; B ennert et a l. 2011) and local inactive galaxies (black triangles; Bennert et al.l 2011), with z = relation (see N2.5|) . The 
errors for the local samples are omitted for clarity (0.4 dex in Mbh, 0.25 dex in Mgph,^). Right panel: The same as in the left panel, 
for total host-galaxy stellar mass. Here, we overplot the 89 BLAGNs from iMerloni et al.l (|201Q ) (blue filled squares; 10 with upper limits 
indicated by arrows). 
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Fig. 3. — Left panel: Offset in log Mbh as a function of constant Mgp^^^ (our objects: red filled pentagons) with respect to the fiducial 
local relation of AGNs (black filled circles). The offset in log Mbh as a function of constant stellar spheroid luminosity from Benn ert et al.l 
(|2010l ) is overplotted (green open symbols), corresponding to AGNs at different redshifts, (left to right: z~0.08, reverberation-mapped 
AGN from Bennert et al. 2010; Bentz et al. 2009; z~0.4 from Bennert et al. 2010; Treu et al. 2007 z~0.6 from Bennert et al. 2010, z^l.8 
from "Bennert et al...201Q : .Peng et al...200& V The best hnear fit derived here is overplotted as dotted fine (Men/Msph,^ oc (1 + z)^-^^^^-^^; 
dashed lines: la ran ^e). Rig;ht panel: T he same as i n the le f t pan el as an offset in log Mbh as a function of constant total host galaxy 
mass (luminosity for lBennert et all 120101 ). The iMerloni et al.l (|2010l ) sample is overplotted (blue filled squares). The lines correspond to 

MBH/Mhost,. OC (1 + Z)1-15±0.15, 
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